Draft version February 21, 2012 

Preprint typeset using fflgX style emulateapj v. 5/2/1 1 



THE 2010 VERY HIGH ENERGY y-RAY FLARE & 10 YEARS OF MULTI-WAVELENGTH OBSERVATIONS OF M 87 

The H.E.S.S. Collaboration: A. Abramowski 1 , F. Acero 2 , F. Aharonian 3,4 ' 5 , A.G. Akhperjanian 6,5 , G. Anton 7 , A. Balzer 7 , A. Barnacka 8 - 9 , 
U. Barres de Almeida 10 , Y. Becherini 11 -' 2 , J. Becker 13 , B. Behera 14 , K. Bernlohr 315 , E. Birsin 15 , J. Biteau 12 , A. Bochow 3 , C. Boisson 16 , 
J. Bolmont 17 , P. Bordas 18 , J. Brucker 7 , F. Brun 12 , P. Brun 9 , T. Bulik 19 , I. Busching 20 - 13 , S. Carrigan 3 , S. Casanova 13 , M. Cerruti 16 , 
P.M. Chadwick 10 , A. Charbonnier 17 , R.C.G. Chaves 3 , A. Cheesebrough 10 , A.C. Clapson 3 , G. Coignet 21 , G. Cologna 14 , J. Conrad 22 , 
M. Dalton 15 , M.K. Daniel 10 , I.D. Davids 23 , B. Degrange 12 , C. Deil 3 , H.J. Dickinson 22 , A. Djannati- Ata'i 1 1 , W. Domainko 3 , L.O'C. Drury 4 , 
G. Dubus 24 , K. Dutson 25 , J. Dyks 8 , M. Dyrda 26 , K. Egberts 27 , P. Eger 7 , P. Espigat 11 , L. Fallon 4 , C. Farmer 2 , S. Fegan 12 , F. Feinstein 2 , 
M.V. Fernandes 1 , A. Fiasson 21 , G. Fontaine 12 , A. Forster 3 , M. Fussling 15 , Y.A. Gallant 2 , H. Gast 3 , L. Gerard 11 , D. Gerbig 13 , B. Giebels 12 , 
JF. Glicenstein 9 , B. Gluck 7 , P. Goret 9 , D. Goring 7 , S. Haffner 7 , J.D. Hague 3 , D. Hampf 1 , M. Hauser 14 , S. Heinz 7 , G. Heinzelmann 1 , 
G. Henri 24 , G. Hermann 3 , J.A. Hinton 25 , A. Hoffmann 18 , W. Hofmann 3 , P. Hofverberg 3 , M. Holler 7 , D. Horns 1 , A. Jacholkowska 17 , 

O.C. de Jager 20 , C. Jahn 7 , M. Jamrozy 28 , I. Jung 7 , MA. Kastendieck 1 , K. Katarzynski 29 , U. Katz 7 , S. Kaufmann 14 , D. Keogh 10 , 
D. Khangulyan 3 , B. Khelifi 12 , D. Klochkov 18 , W. Kluzniak 8 , T. Kneiske 1 , Nu. Komin 21 , K. Kosack 9 , R. Kossakowski 21 , H. Laffon 12 , 
G. Lamanna 21 , D. Lennarz 3 , T. Lohse 15 , A. Lopatin 7 , C.-C. Lu 3 , V. Marandon 11 , A. Marcowith 2 , J. Masbou 21 , D. Maurin 17 , N. Maxted 30 , 
M. Mayer 7 , T.J.L. McComb 10 , M.C. Medina 9 , J. Mehault 2 , R. Moderski 8 , E. Moulin 9 , C.L. Naumann 17 , M. Naumann-Godo 9 , M. de Naurois 12 
D. Nedbal 31 , D. Nekrassov 3 , N. Nguyen 1 , B. Nicholas 30 , J. Niemiec 26 , S.J. Nolan 10 , S. Ohm 32 - 253 , E. de Ona Wilhelmi 3 , B. Opitz 1 , 
M. Ostrowski 28 , I. Oya 15 , M. Panter 3 , M. Paz Arribas 15 , G. Pedaletti 14 , G. Pelletier 24 , P.-O. Petrucci 24 , S. Pita 11 , G. Puhlhofer 18 , 
M. Punch 11 , A. Quirrenbach 14 , M. Raue 1 -*, S.M. Rayner 10 , A. Reimer 27 , O. Reimer 27 , M. Renaud 2 , R. de los Reyes 3 , F. Rieger 3,33 , J. Ripken 22 
L. Rob 31 , S. Rosier-Lees 21 , G. Rowell 30 , B. Rudak 8 , C.B. Rulten 10 , J. Ruppel 13 , V. Sahakian 6,5 , D.A. Sanchez 3 , A. Santangelo 18 , 
R. Schlickeiser 13 , F.M. Schock 7 , A. Schulz 7 , U. Schwanke 15 , S. Schwarzburg 18 , S. Schwemmer 14 , F. Sheidaei 11,20 , J.L. Skilton 3 , H. Sol 16 , 

G. Spengler 15 , L. Stawarz 28 -*, R. Steenkamp 23 , C. Stegmann 7 , F. Stinzing 7 , K. Stycz 7 , I. Sushch 15 , A. Szostek 28 , J.-P. Tavernet 17 , 
R. Terrier 11 , M. Tluczykont 1 , K. Valerius 7 , C. van Eldik 3 , G. Vasileiadis 2 , C. Venter 20 , J.P Vialle 21 , A. Viana 9 , P. Vincent 17 , H.J. Volk 3 , 
F. Volpe 3 , S. Vorobiov 2 , M. Vorster 20 , S.J. Wagner 14 , M. Ward 10 , R. White 25 , A. Wierzcholska 28 , M. Zacharias 13 , A. Zajczyk 82 , 
AA. Zdziarski 8 , A. Zech 16 , H.-S. Zechlin 1 , 

AND 

The MAGIC Collaboration: J. Aleksic 34 , L. A. Antonelli 35 , P. Antoranz 36 , M. Backes 37 , J. A. Barrio 38 , D. Bastieri 39 , J. Becerra 
Gonzalez 40 , W. Bednarek 41 , A. Berdyugin 42 , K. Berger 40 , E. Bernardini 43 , A. Biland 44 , O. Blanch 34 , R. K. Bock 45 , A. Boller 44 , 
G. Bonnoli 35 , D. Borla Tridon 45 , I. Braun 44 , T. Bretz 46 , A. Canellas 47 , E. Carmona 45 , A. Carosi 35 , P. Colin 45 *, E. Colombo 40 , 
J. L. Contreras 38 , J. Cortina 34 , L. Cossio 48 , S. Covino 35 , F. Dazzi 48 , A. De Angelis 48 , E. De Cea del Pozo 49 , B. De Lotto 48 , C. Delgado 
Mendez 40 , A. Diago Ortega 40 , M. Doert 37 , A. Dominguez 50 , D. Dominis Prester 51 , D. Dorner 44 , M. Doro 52 , D. Elsaesser 46 , D. Ferenc 51 , 
M. V. Fonseca 38 , L. Font 52 , C. Fruck 45 , R. J. Garcia Lopez 40 , M. Garczarczyk 40 , D. Garrido 52 , G. Giavitto 34 , N. Godinovic 51 , D. Hadasch 49 
D. Hafner 45 , A. Herrero 40 , D. Hildebrand 44 , D. Hohne-Monch 46 , J. Hose 45 , D. Hrupec 51 , B. Huber 44 , T. Jogler 45 , S. Klepser 34 , 
T. Krahenbuhl 44 , J. Krause 45 , A. La Barbera 35 , D. Lelas 51 , E. Leonardo 36 , E. Lindfors 42 , S. Lombardi 39 , M. Lopez 38 , E. Lorenz 44 , 
M. Makariev 53 , G. Maneva 53 , N. Mankuzhiyil 48 , K. Mannheim 46 , L. Maraschi 35 , M. Mariotti 39 , M. Martinez 34 , D. Mazin 34 -*, M. Meucci 36 , 
J. M. Miranda 36 , R. Mirzoyan 45 , H. Miyamoto 45 , J. Moldon 47 , A. Moralejo 34 , P. Munar 47 , D. Nieto 38 , K. Nilsson 42 , R. Orito 45 , I. Oya 38 , 
D. Paneque 45 , R. Paoletti 36 , S. Pardo 38 , J. M. Paredes 47 , S. Partem 36 , M. Pasanen 42 , F. Pauss 44 , M. A. Perez-Torres 34 , M. Persic 48 , 
L. Peruzzo 39 , M. Pilia 54 , J. Pochon 40 , F. Prada 50 , P. G. Prada Moroni 55 , E. Prandini 39 , I. Puljak 51 , I. Reichardt 34 , R. Reinthal 42 , W. Rhode-" 
M. Rmo 47 , J. Rico 56 , S. Rugamer 46 , A. Saggion 39 , K. Saito 45 , T. Y. Saito 45 , M. Salvati 35 , K. Satalecka 43 , V. Scalzotto 39 , V. Scapin 38 , 
C. Schultz 39 , T. Schweizer 45 , M. Shayduk 45 , S. N. Shore 55 , A. Sillanpaa 42 , J. Sitarek 41 , D. Sobczynska 41 , F. Spanier 46 , S. Spiro 35 , 
A. Stamerra 36 , B. Steinke 45 , J. Storz 46 , N. Strah 37 , T. Suric 51 , L. Takalo 42 , H. Takami 45 , F. Tavecchio 35 , P. Temnikov 53 , T. Terzic 51 , 
D. Tescaro 55 , M. Teshima 45 , M. Thom 37 , O. Tibolla 46 , D. F. Torres 56 , A. Treves 54 , H. Vankov 53 , P. Vogler 44 , R. M. Wagner 45 , Q. Weitzel 44 , 

V. Zabalza 47 , F. Zandanel 50 , R. Zanin 34 , 

AND 

The VERITAS Collaboration: T. Arlen 57 , T. Aune 58 , M. Beilicke 59 *, W. Benbow 60 , A. Bouvier 58 , S. M. Bradbury 61 , J. H. Buckley 59 , 
V. Bugaev 59 , K. Byrum 62 , A. Cannon 63 , A. Cesarini 64 , L. Ciupik 65 , M. P. Connolly 64 , W. Cui 66 , R. Dickherber 59 , C. Duke 67 , M. Errando 68 , 
A. Falcone 69 , J. P. Finley 66 , G. Finnegan 70 , L. Fortson 71 , A. Furniss 58 , N. Galante 60 , D. Gall 72 , S. Godambe 70 , S. Griffin 73 , J. Grube 65 , 
G. Gyuk 65 , D. Hanna 73 , J. Holder 74 , H. Huan 75 , C. M. Hui 70 *, P. Kaaret 72 , N. Karlsson 71 , M. Kertzman 76 , Y. Khassen 63 , D. Kieda 70 , 

H. Krawczynski 59 , F. Krennrich 77 , M. J. Lang 64 , S. LeBohec 70 , G. Maier 78 , S. McArthur 59 , A. McCann 73 , P. Moriarty 79 , R. Mukherjee 68 , 
P. D Nunez 70 , R. A. Ong 57 , M. Orr 77 , A. N. Otte 58 , N. Park 75 , J. S. Perkins 80 81 , A. Pichel 82 , M. Pohl 83 78 , H. Prokoph 78 , K. Ragan 73 , 

L. C. Reyes 75 , P. T. Reynolds 84 , E. Roache 60 , H. J. Rose 61 , J. Ruppel 83 - 78 , M. Schroedter 60 , G. H. Sembroski 66 , G. D. §enturk 68 , 

I. Telezhinsky 83 ' 78 , G. Tesic 73 , M. Theiling 60 , S. Thibadeau 59 , A. Varlotta 66 , V. V. Vassiliev 57 , M. Vivier 74 , S. P. Wakely 75 , T. C. Weekes 60 , 

D. A. Williams 58 , B. Zitzer 66 , 

AND 

U. Barres de Almeida 45 , M. Cara 85 , C. Casadio 86 - 50 , C.C. Cheung 87 , W. McConville 88 - 89 , F. Davies 90 , A. Doi 91 , G. Giovannim 86 - 92 , 
M. Giroletti 86 , K. Hada 93 - 94 , P. Hardee 95 , D. E. Harris 96 , W. Junor 97 , M. Kino 94 , N.P. Lee 96 , C. Ly 98 , J. Madrid 99 , F. Massaro 96 , 
C. G. Mundell 100 , H. Nagai 94 , E. S. Perlman 85 , I. A. Steele 100 , R.C. Walker 101 , D.L. Wood 102 

Draft version February 21, 2012 

ABSTRACT 

The giant radio galaxy M87 with its proximity (16Mpc), famous jet, and very massive black hole ((3 - 
6) x 10 M ) provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) 
y-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has 
been established as a VHE y-ray emitter since 2006. The VHE y-ray emission displays strong variability on 
timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the 
MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected 
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triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent 
sampling of the VHE y-ray light curve enables one to derive a precise temporal characterization of the flare: the 
single, isolated flare is well described by a two-sided exponential function with significantly different flux rise 
and decay times of T™ e = (1.69 + 0.30) days and T ^ cdy = (0.611 ± 0.080) days, respectively. While the overall 
variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 
2008, they share very similar timescales (~day), peak fluxes (fl>>o.35TeV - (1 — 3) X 10~ H phcirT 2 s _1 ), and 
VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in 
contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided 
with a VHE flare. On the other hand, Chandra X-ray observations taken ~ 3 days after the peak of the VHE 
y-ray emission reveal an enhanced flux from the core (flux increased by factor ~ 2; variability timescale < 2 
days). The long-term (2001-2010) multi-wavelength (MWL) light curve of M 87, spanning from radio to VHE 
and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE y-ray 
emission. No unique, common MWL signature of the three VHE flares has been identified. In the outer kpc 
jet region, in particular in HST-1, no enhanced MWL activity was detected in 2008 and 2010, disfavoring it 
as the origin of the VHE flares during these years. Shortly after two of the three flares (2008 and 2010), the 
X-ray core was observed to be at a higher flux level than its characteristic range (determined from more than 
60 monitoring observations: 2002-2009). In 2005, the strong flux dominance of HST-1 could have suppressed 
the detection of such a feature. Published models for VHE y-ray emission from M 87 are reviewed in the light 
of the new data. 

Subject headings: galaxies: active - galaxies: individual (M 87) - gamma rays: galaxies - galaxies :jets; nuclei 
- radiation mechanisms: non-thermal 
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1. INTRODUCTION 

The giant radio galaxy M 87 provides a unique environment 
to study relativistic plasma outflows and the surroundings of 



cussed as a possible VHE y-ray emission site (e.g. Stawarz 



su per-m assive black holes (SMBH). Its prominent jet (Cur 
|tis||1918 1 is resolved from radio to X-r ays displaying com 



plex st ructures (knots, diff use emission; |Perlman et al.|1999 



|2001a) >, strong variability ^Harris et al.|2003| |2006|l, and ap 
parent super-luminal motion (Biretta et al. 1999 1 Cheung et al. 
[2007) . With its proximity (16.7 + 0.2 Mpc; |Mei et al 2007 ) 
and its very massive black hole of M B h - (3 - 6) X 10 9 M e ™ 
acchetto et al.|1997| IGebhardt & Thomas|200"9[ ) high res- 
ton very long baseline interferometry (VLBI) at radio 
wavelengths enables one to directly probe structures with 
sizes down to < 200 Schwarzschild radii. From the detection 
of super-luminal features in the jet in the optical and radio the 
jet orientation angle towards t he line of sight at the su b-kpc 
scale is limited to 6 < 20° ( |Biretta et al.|1999[ [Cheung et ah 
[2007). 

Evidence for very high energy (VHE; E > 100 GeV) y-ray 
emission from M 87 was reported by the HEGRA collabo 
ration in 2003 ( Aharonian et al.|[200 3 | 105 | an d was later con 



firmed by H.E.S.S., VERITAS, and MAGI C (Aharonia n et al. 
|2006a| lAcciari et al.||2008] |Albert et al.||2008a| ). While the 



large majority of active galactic nuclei (AGN) with detected 
VHE y-ray emission are strongly beamed sources, M 87 is one 
of only four known VHE AGN with weak or, at most, mod- 
erate bea ming; the other three being the radio galaxies Cen - 
taurus A ( jAharonian et al.|2009|l, IC 310 flAleksiB et al.|20"Tol >, 
and NGC 1275 ( |Mariotti & MAGIC Collaboration^ 10a| l. In- 
terestingly, for such a weakly beamed source, M87 shows 
strong variability at VHE with timescales of the order of days 
(Ahar onian et al.||2006a| |Albert et "aL"1|2008a| |Acciari et al.| 
|2009) >. This points, through the causality argument, towards a 
compact emission region < 5 x 10 15 5cm (6 being the Doppler 
factor of the emitting plasma) corresponding to only a few 
Schwarzschild radii R s = 2GM Bii /c 2 10 15 cm. At GeV en- 
ergies M 87 has recently been detected as a weak source by 
Fenni-LAT ( |Abdo et al.|2009] >. 

The exact location of the VHE y-ray emitting region in 
M 87 remains elusive. The angular resolution of ground-based 
VHE instruments is of the order of 0.1° (corresponding to 
~ 30kpc projected size) and, therefore, does not allow for a 
direct precise determination of the VHE y-ray emission site in 
the inner kpc-scale structures, alth ough the outer radio lobes 
can be excluded as the origin (Aharonian et al. 2006a| l. To 
further investigate the location of the VHE y-ray emission site 
and the associated production mechanisms, variability studies 
and the search for correlations with other wavele ngths have 
successfully been utilized (e.g. |Acciari et al.|2009] ). Of par- 
ticular interest are radio observations, since they allow for 
the highest angular resolution, and X-ray observations, due 
to their potential connection with the VHE y-ray emission in 
e.g. synchrotron self-Compton (SSC) models. 

Up to now, three episodes of enhanced VHE y-ray emis- 
sion have been detected from M 87, with details on the lat- 
est one, observed in 2010, being reported in this paper. The 
first one, detected in 2005 ( |Aharonian et al.|20 06a), coincided 
with a n extreme multi-frequency outburst of the jet feature 
HST-1 dHarris et al.^OOSlfeOOelF 5 ! which has also been dis- 



In the followi ng Mrh - 3 X 10 9 M Q is adopted 



105 See also 



Le Bohec et al. 



(2004 



et al 
the 



1. 2006| |Cheung et al |[2007] fHarris et al ||2009 
second flaring episode, detected in 2008, HST 



Durins 



HST-1 was in 

a low flux state, but radio measurements at 43 GHz with the 
Very Long Baseline Array (VLBA) showed a flux increase in 
the core region within a few hundred Schwarzschild radii of 
the SMBH, suggesting the direct vi cinity of the SM BH as the 
origin of the VHE y-ray emission (Acc iari et al.|20 09). This 
conclusion was further supported by the detection of an en- 
hanced X-ray flux from the core region by Chandra. 

In this paper results from a joint observation campaign of 
M 87 in 2010 are presented and discussed in the broader con- 
text of the multi-wavelength (MWL) behavior of M 87 over 
the past ten years. During the campaign, a high flux state at 
VHE was detected (|Mariotti & MAGIC Collaboration^ 10b| 
|Ong & Mariotti||2010b . Characteristics of the VHE flare are 
investigated and possible correlations with other wavelengths 
are discussed. New observational results from H.E.S.S., 
MAGIC, VERITAS, Fenni-LAT, Chandra, HST, LT, VLBA, 
MOJAVE, VLA, and the EVN are presented^ 

In Sec. [2] the instruments and the data are introduced. In 
Sec. [3] the characteristics of the VHE high state are investi- 
gatedand compared to previous flares. New results from opti- 
cal polarimetry observations with the Hubble Space Telescope 
and the Liverpool Telescope are presented and discussed in 
Sec. |4] In Sec. [5] the VHE results are confronted with the 
broader MWL picture and the theoretical implications con- 
cerning the VHE y-ray emission site are presented in Sec. [6] 
The paper concludes with Sec. [7] 

2. DATA 

2.1. Very high energy (VHE) 

M 87 has been the target of several coordinated VHE mon- 
itoring campaigns since 2008. In 2010 the source was jointly 
monitored by MAGIC and VERITAS. Duri ng the campaign 
Mariotti & MAGIC Colla boration (2010b) reported an in- 
creased flux from the source in February 2010, but follow-up 
observations did not reveal further activity at VHE. A sec- 
ond increase in flux, marking the onset of a strong VH E flare, 
was detected in April 2010 by MAGIC and VERITAS (|Ong & 



Mario tti|2010| ) triggering further ToO (Target of Opportunity) 
observations by H.E.S.S., Chandra, the VLBA, and the EVN. 

For the 2010 campaign, M 87 has been observed for a total 
of ~80h from December 2009 to June 2010. The data from 
different VHE instruments have been combined after separate 
analysis within the individual collaborations. Integral fluxes 
for the VHE band <J>vhe are calculated above an energy of 
350 GeV. Observations taken or published with a different en- 
ergy threshold are extrapolated using the average measured 
VHE spectrum, which is well described by a power la w spec- 
trum djV/dE oc £- r with photon index r = 2.3 ± 0.11 ( |Albert| 
|et al.|2008a| ). 

Given the indications for spectral variability of M87 at 



s pe 

VHE JAharonian et al.|2006al [Albert et al.|2008aHAliu et al 



201 1\ using a single photon index for the flux extrapolation 
could introduce a bias in the light curves. For the flaring 
states the hard est spectral index reported is F = 2.21 + 0.18 
( | Albert et~a l. 2008a) and, therefore, the systematic error in- 
troduced when using F = 2.3 is completely negligible com- 
pared to the typical statistical error of 10-20% of the flux dur- 
ing such states. For the quiescent state spectral indices up 



The outburst was also followed by the ejection of apparent superlumi- 



nal radio components from HST-1 lCheung ~et al.| 

107 Full names of the instruments can be found in Sec. [5] 
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Fig. 1. — Multi-wavelength light curve of M 87 from 2001 to 201 1. The VHE y-ray flux (top panel) is calculated above an energy threshold of 350 GeV (see 
text). Separate fluxes for the core and HST-1 are shown in cases where the instrument resolution is sufficient to separate the two components. Gray vertical bands 
mark the times of increased VHE activity in 2005, 2008, and 2010 (see Fig. pi. The dashed line and the gray horizontal band in the 2nd panel marks the average 
flux with 1 s.d. error measured by Fermi-LJXT, The radio flux of HST-1 at (Efferent frequencies has been normalized to the 5 GHz flux assuming a spectrum 
S v ~ v~" with a = 0.6. All flux errors shown are the 1 s.d. statistical errors except for the LT data where the uncertainty on the contribution from the galaxy is 
included in the error bars. For details on the data, the data analysis, and references see text. 



to T = 2.60 ± 0.30 have been reported (Albert et al. 2008a). 
This could, in principle, create a bias for the H.E.S.S. data set, 
where the flux is extrapolated down to 350 GeV from energies 
above 500 GeV. For example, for spectral indices of F t = 2.3 
and F2 = 2.6 the flux extrapolated down to 350 GeV from 
700 GeV would differ by a factor (350 GeV/700 GeV) r2 ~ ri ~ 
0.8, which would imply errors of order ~20%. Given that 
the typical flux error for the nightly averaged flux bins in the 
quiescent state is of order ~100% (i.e. the individual nightly 
flux points are not significant detections) the systematic un- 
certainty introduced by using a single photon index for ex- 
trapolation can also safely be neglected. 



Given recent results on the Crab Nebula (the reference 
source for ground-based VHE instruments) indicating only 
small systematic offset between the en ergy scale of differ- 
ent instruments (M eyer et al.| [2010 108 and the general good 
agreement of the flux measurements during quasi simultane- 
ous observations, the systematic error between the different 
instruments is estimated to be small compared to the statisti- 
cal error of individual measurements of the 2010 campaign. 

Additional archival data from |Aharonian et al. | (f2006a)> , | Ac- 1 
|ciari et al.|p008) , [Albert et al.| ( |2008a| l, |Acciari et al.| ( |2010| ), 

108 Energy scale scaling factors relative to the Fermi-LAT energy scale of 
0.961 ± 0.004 for H.E.S.S. and 1.03 ± 0.01 for MAGIC are derived. 
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and |Aleksic et al. (201 la I are also shown in the light curves 
(Fig.mandFig.|2f 

In the following, the main characteristics of the VHE ob- 
servatories involved in the M 87 campaign are reviewed, and 
details on the corresponding data-sets are presented. 



indicating a spe ctral change w ith flux will be published in a 
parallel paper ( Aliu et al.|201 1 1. 



H.E.S.S, 
(H.E.S.S 



109 



- The High Energy Stereoscopic System 
is a system of four large (13 m mirror di- 
ameter) imaging atmospheric Cherenkov telescopes (IACTs) 
for the detection of VHE y-rays, located in the southern 
hemisphere in Namibia (23°16'S, 16°30'W; 1800 m above 
sea level). It has been in operation sinc e 2002, with the full 
array completed in 2004 ( |Hinton|[2004| . H.E.S.S. observed 
M87 for 10.7 h in 2010 (dead-time corrected) yielding a 
total detection signifi cance of 9.7 standard deviations (s.d.; 
following |Li & Ma| |1983|l using a standard H illas-type 
analysis and cuts fromlAharonian et al. 2006b| (software 
version hap- 18-06). M 87 culminates at ~35° zenith angle 
at the H.E.S.S. site resulting in an energy threshold of 
Zithr > 500 GeV for the analysis and data-set. The integral 
flux is extrapolated down to 350 GeV using the average 
energy spectrum (see above). In addition, results from a 
re-analysis of the 2004 and 2006 H.E.S.S. data in nightly 
flux bins, utilizing the same analysis as discussed above, 
are presented in the light curves. Cross check on the results 
with data from an independent calibration chain and utilizing 
different analysis have been performed and good agreement 
is found. 

MAGIC — The M ajor Atmospheric Gamma-ray Imaging 



Cherenkov (MAGIC ) 110 | telescope system consists of two 
17 m diameter IACTs located at the Roque de los Muchachos 
Observatory, on the Canary Island of La Palma (28°46' N, 
17°53' W; 2200 m above sea level). Since 2005 M 87 h as been 
regularly observed by MAGIC with a single telescope ( Albert 
|etaLl|2008b[ |Aleksic et aEl|2011a) . In the Autumn of 2009 
the system became stereoscopic, with the commissioning of 
a second telescope and the stereo trigger, resulting in an al- 



most doubling of its sensitivity (|Colin et al. 2009 Aleksic 



|et al.||201 lb] ). In 2010, M87 observations were conducted 
for the first time, in stereoscopic mode. 20 h of good quality 
data were taken between January and June with zenith angles 
ranging from 16 to 35 degrees . Analysis of these d ata with 
the standard MAGIC software (Moralej cTet al.|2009| resulted 
in a 10 s.d. detection above 200 GeV 

VERITAS — The Very Energeti c R adiation Imaging Tele- 
scope Array System (VERITAS^ consists of four 12 m 
diameter IACTs and is located at the base camp of the 
Fred Lawrence Whipple Observatory in southern Arizona 
(31°40'N, 110°57' W; 1280m above sea level). More details 
about VERITAS, the data calibration, and th e analysis tech- 
niques can be found in |Acciari et al. ( 2008) 1. VERITAS ob- 
served M87 in 2010 for 48.2 h (after quality cuts) resulting 
in a detection of 26 s.d. The zenith angles of the observations 
ranged from 19 to 40°, with a few nights with zenith angles up 
to 60° during the flares (April 9-1 1th). The energy threshold 
of the analysis for the mean zenith angle of the observations of 
25° is 250 GeV. The VERITAS data covered the whole 2010 
flare period; a detailed study of the VHE spectral evolution 



2.2. High energy (HE) 

Fermi-IAT — The Large Area Telescope (LAT) onboard 
the Fermi satellite is a pair-conversion telescope that covers 
the e nergy range from 2 MeV to more than 300 GeV 
(HE) ( |Atwood et"aLl|2009"l >. The LAT instrument features a 
per-photon angular resolution of 6* 6 g% = 0.8° at 1 GeV and a 
large field-of-view of 2.4 sr. The primary mode of operation 
is an all-sky survey mode, where the full sky is covered ap- 
proximately every three hours. The LAT data for this analysis 
consists of two years of nominal all-sky survey data between 
the energy range 100 MeV and 300 GeV, and spanning the 
mission elapsed time (MET) 239557417 to 302630530 
(August 4, 2008 through August 4, 2010). Event selections 
include "diffuse" class events recommended for point source 
analysis, a rocking angle cut of < 52°, and a zenith angle cut 
of < 100° in order to avoid contamination from the Earth's 
limb. A 2 ks window beginning at MET 259459364 was also 
removed in order to avoid contamination from GRB 090323, 
which occurred nearby. All LAT analysis was performed 
using instrument response functions (IRFs) P6_V11_DIFFUSE 
and science tools v9r2Qp0, along with the recommendec 112 
Galactic diffuse gll_iem_v02_P6_VHJ)IFFUSE . fit 
and corresponding isotropic spectral template 
isotropic J. em_v02_P6_Vll_DIFFUSE.txt. 

An analysis of the LAT spectrum over the two-year period 
was perform ed using a binned likelihood analysis (Mattox 
et al.||1996] > selecting all events that fell within a 20 x 20 1 



109 http://www.mpi-hd.mpg.de/hfm/HESS/ 

1 10 http://wwwmagic.mppmu.mpg.de/ 

1 1 1 http://veritas.sao.arizona.edu/ 



square region of interest (ROI) centered at the M 87 radio po- 
sition of the core. All point sources from an internal two-year 
preliminary catalog that fell within 15° of the source were 
included in the fit. All sources within the square ROI were 
modeled with a power law spectrum with their normalization 
and index as free parameters, while those that fell outside of 
the ROI were fixed to their catalog values. The M 87 spectrum 
was modeled as a power law with photon index and normal- 
ization para meters left free and using the radio position ( |Fey| 
|et al.|[2 004) as the source locati on. A point source was de- 
tected with a test statistic (T S; |Mattox et al.|[T9"96] l of 301, 
representing a detection of V301 ^ 17 s.d. From the result- 
ing fit, the photon index and flux (> 100 MeV) were found 
to be 2.16 ± 0.07 and (2.66 ± 0.36) x lO^phcm^s 1 , re- 
spectively. In order to test for curvature, the spectrum was 
also fit to a log parabola, where the TS of the overall fit was 
found to improve by only 0.89, which does not represent a 
statistically significant improvement over the single power 
law. The largest systematic errors can be attributed to un- 
certainties in the modeling of the diffuse background emis- 
sion. These were estimated by repeating the analysis with 
both binned and unbinned gtlike using a refined version 
of the current diffuse background model that is under devel- 
opment by the LAT collaboration. Systematic errors on the 
index and flux were thus found to be (+0.05/ - 0.01) and 
(+0.40/ - 0.13) x 10 _8 phcm _2 s _1 , respectively. Compar- 
ing results from this analysis for the last 14 month with the 
initial 10-month spectrum reported in |Abdo et al.| ( [2009) l, no 
evidence of variability in the flux above 100 MeV was found. 
Comparing the flux above 1 GeV between these two epochs, 
however, a marginal indication of a rise in the flux of the latter 

112 http : //f ermi . gs£c . nasa . gov/ssc/data/access/lat/ 

Backgr oundMode 1 s . html 



2010 VHE flare & 10 years of MWL observations of M 87 



7 



epoch at a significance of 2s.d. is found. Systematic uncer- 
tainties on the LAT flux due to the instrument response are 
estimated at values of 10 %, 5 %, and 20 % above and below 
their n ominal values at log (E/MeV) = 2, 2.75, and 4, respec- 
tively ( Atwoo d et al.|2009| > . 

The LAT light curve (> 1 GeV) was constructed using the 
events that fell within a 10° circular ROI centered at the M 87 
radio position. Fluxes were derived for bins with a detec- 
tion significance exceeding 3 s.d., otherwise upper limits were 
calculated. To generate the light curve, events were grouped 
into 56-day time bins, and a separate likelihood analysis us- 
ing gtlike was performed over each of the bins. All point 
sources from the two-year fit were included in the model over 
each interval. Sources that fell within the 10° ROI were fit 
with their normalization parameters free, while the photon in- 
dex of each source was fixed to the best-fit value obtained 
from the full two-year analysis. Both the index and normal- 
ization parameters of M 87 were left free, except in the case 
of upper limit calculations, in which case the photon index 
was fixed to the nominal two-year average of 2.16. In order 
to avoid modeling sources with a negative TS value, an initial 
fit over each interval was performed, and all sources found to 
have a TS< 1 were subsequently removed from the fit . Fol- 
lowing the method for variability detection outlined in |Abdo| 
et al. (2010a), the weighted average flux was first calculated 



with a resulting value of ( 1 .62 ± 0. 1 8) x 1 0~ 9 ph cirT 2 s~ 1 . A x 2 
analysis was then performed by comparing the best-fit values 
of all points to the weighted average, and the resulting prob- 
ability P(x 2 > xlbJ was f° unc l t0 be 0.027, which represents 
a significance of 2.2 s.d. for the source to be variable an d falls 
slightly below the threshold for variability defined in Abdo 
[etaLl ( |2010a| l. 



2.3. X-ray 

Chandra — X-ray data have been taken with the Advanced 
CCD Imaging Spectrometer (ACIS) on board the Chandra 
satellite. For details of the Chandra data reduction proce- 
dures, see |Harris et all ((2003), |Perlman et aE] ( |2003| ), and 



Harris et al.| ( 2009 1. In brief, a l/8th segment of the back 
illuminated S3 chip of the ACIS detector aboard Chandra is 
used. This results in a frame time of 0.4s with 90 % effici ency. 
Although thi s setup was essentially free of pileup when Wil- 
son & Yang ( 2002 ) tested various options during 2000 July, 
with the advent of the ever increasing brightness of HST-1, 
pileup ( Davis|2001[ l became a major problem so the measure 
of intensity was switched to a detector based unit: keV/s. This 
approach uses the event 1 file with no grade filtering (so as to 
recover all events affected by 'grade migration') and energies 
from 0.2 to 17 keV are integrated so as to recover all the en- 
ergy of the piled events. Other uncertainties for piled events 
come from the on-board filtering, the 'eat-thy-neighbor' ef- 
fect, and second order eff ects such as rele ase of trapped charge 
(see Appendix A of Harris et al. 2009[ ). Although a small 
circular apertur e was used for flux-map photometry in Har- 
ris et al.| ( |2003j ), the basic analysis for this pap er adopts the 



rectangular regions used in Harris et al. (2006) so as to en 
compass more of the point spread function (PSF). All events 
within each rectangle are weighted by their energy and the 
sum of these energies, when divided by the exposure times, 
gives the final keV/s value used in the light curve. Uncer- 
tainties are strictly statistical, based on the number of counts 
measured ( "VTY /AO and typically range from 1% to 5%. 



2.4. Optical 

Hubble Space Telescope (HST) — Near Ultraviolet HST ob- 
servations were obtained with the Space Telescope Imaging 
Spectrograph (STIS) and the Advanced Camera for Surveys 
High Resolution Camera (ACS/HRC). STIS imaging was 
taken with the F25QTZ filt er that has its maximum through- 
put wavelength at 2364.8 A ( |Kim Quijano et al.|2003| . All ob- 
servations after August 2004 were taken with the ACS/HRC 
using two filters: F220W and F250W which have their max- 
imum throughp ut wavelength at 225 5.5 A and 27 15.9 A, re- 
spectively (Mack & Gilliland 2003). All science ready files 
were retrieved from the STScl public archive and processed 
throug h the PYRAF task MULTIDRIZZEL ( |Fruchter etak 
2009 ). Both HST detectors have a pixel scale of ~ 0.024" per 
pixel and a resolution that enables one to clearly separate the 
nucleus and the innermost components of the jet, i.e. H ST- 1 . 
The details of these observations have been presented in|Perl- 
man et al.| ( |2003l ), |Madrid| ( |2009l ), and |Perlman et al.| ( |20TT] L 

Optical HST polarimetry observations were obtained with 
the ACS/HRC and the Wide-field Planetary Camera 2 
(WFPC2), using the F606W filter that has its maximum 
throughput at roughly ~6000 A and a nearly flat throughput 
curve from 4800 to 7000 A. The single WFPC2 observation 
(pixel scale ~ 0.1 arcsec/pixel) was obtained when the ACS 
was not operational, and used the POLQ polarizing filter, 
while the ACS/HRC observations used the POLVIS polariz- 
ing filter. For both of these polarizers, the F606W filter gives 
nearly optimal transmission of parallel polarized light as well 
as rejection of perpendicularly polarized light. Polarimetry 
observations were obtained at 18 epochs between December 
2002 and November 2007, with 14 of the 18 observations con- 
centrated between November 2004 and December 2006 on 
the same schedule as the ultraviolet photometry. As with the 
UV observations, all polarimetry was obtained from the HST 
archive and re-calibrated using updated flat field files and im- 
age distortion correct ion (IDC) tables (Mob asher et al.|2002 



Pavlovsky et al. 2002 \ . multidrizzle was used to combine 



and cosmic- ray reject the images, which were aligned using 
Tweakshifts (Frucht er et al.|2009) , and CTE (Charge Transfer 
Efficiency) corrections were computed using d ata found in the 
ACS instrument handbook (Boffi et al. 2007 ). Once this was 
done, the images in each polarizer were combined according 
to recipes i n the ACS and WFPC2 instrument han dbooks, re- 
spectively ( [Boffi et al.|2007)|Biretta & McMaster| 1997b . Be- 
fore performing photometry and polarimetry, galaxy emission 
was subtracted from the images. This was done using a model 
computed in the Stokes I image using the IRAF tasks EL- 
LIPSE and BMODEL. Error bars for the polarimetry are typ- 
ically 2-3% for high signal-to- noise data (further details can 
be found in Perlma net al.|201l} . 

Liverpool Telescope (LT) — Hybrid R+V-band optical po- 
larimetry data (460-720 nm at F WHM ) were taken with the 
2-meter Liverpool Telescope (LT 113 located on La Palma, 
using the newly commissioned R1N002 fast-readout imaging 
polarimeter ( |Steele et al.[2 010). The polarimeter uses a rotat- 
ing polaroid with frequency of approximately 1 Hz, during the 
cycle of which 8 exposures of the source are obtained. These 
exposures are synchronized with th e phase of the polaroid and 
following the analysis method of Clarke & Neumayer ( 2002[ ) 
allows determination of the degree and angle of polarization. 

113 http://telescope.livjm.ac.uk/ 
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The LT observations were taken during and shortly after 
the observed VHE high-state, MJD 55295-402, with three 
measurements taken contemporaneously to the time of the 
VHE flare, MJD 55295-97. Further data on M 87 have been 
taken since then and will be presented elsewhere. Total in- 
tegration times of typically 100 s were used in the observa- 
tions, corresponding to an achieved polarimetric accuracy of 
about 1 % for the brightness of the source. The field of view 
of the instrument is of 4 x 4arcmin with a pixel scale of 
0.45 arcsec/pixel, which at the distance of M 87 corresponds 
to a linear scale of ~40pc/pixel, equivalent to approximately 
half the distance between the core and the innermost jet com- 
ponent, the knot HST-1. Given the typical seeing during the 
observations of 1.0" to 1.7" (FWHM) the two components, 
core and HST-1, cannot be resolved individually. There- 
fore, an aperture radius for the integration of the signal of 
2.7 arcsec was used, so that both the nucleus and HST-1 were 
included. The outer jet is nevertheless well resolved and in- 
dependent light-curves were produced with the same aperture 
radius used for the core but now centered at knot-A, located 
12.3 arcsec downstream from the nucleus, revealing a steady 
polarization aligned with the jet direction. 

The greatest source of error in the determination of the po- 
larization levels is a systematic effect caused by contamina- 
tion by the bright host galaxy, which extends out to several 
kpc and dilutes the polarization signal from the nucleus. The 
strength of the contamination is estimated by adding a Sersic 
profile with n — 4 to a frame built by summing data from all 
polaroid angles. This sho ws the well known flat tening in the 
central ~ 10 arcsec (e.g. Kormend y et al.|2009| . Fitting the 
data from 9 to 2.7 arcsec reveals an excess within the central 
2.7 arcsec due to the core plus HST-1 of ~ 10 + 5%. The 
photometric flux of the background galaxy is of order 90 % of 
the total flux within the aperture used for extracting the polar- 
ization measurements. The measured fractional polarization 
was therefore multiplied by a factor of 10 and the flux mea- 
surement divided by a factor 10 in order to correct for this 
background light, and allow comparison with the high resolu- 
tion HST measurements of the core + HST-1. 

Light-curves for the polarization position angle (or electric 
vector position angle, EVPA, defined as 0.5 x aictan(u/q), 
where u and q are the linear Stokes parameters) were ob- 
tained for all the epochs of observation, and are presented in 
Fig. [4] along with the other polarization quantities. Due to 
host-galaxy contamination,the polarization position angle for 
the nucleus and HST-1 combined are measured to an accu- 
racy of 25°, considerably degraded in comparison to the 10° 
resolution achieved for the outer jet where the host galaxy is 
fainter. 

2.5. Radio 

Radio interferometers enable one to observe the jet of M 87 
with a large variety of angular resolutions and sensitivities, 
depending on the array size and observing frequency. In 
Tab. [T]the highest angular resolution achieved for each instru- 
ment contributing data to this paper is shown. Overall, they 
span an angular scale range from as large as ~ 0.1" (from 
the VLA at 22 GHz) down to a fraction of a milli-arcsecond 
(e.g. ~ 0.2 mas with the VLBA at 43 GHz, see Table 1). 
The longest baseline and highest frequency observations pro- 
vide the most valuable information about the compact, flat 
spectrum core; conversely, instruments with shorter baselines 
(such as the VLA) or lower observing frequency (like the 
EVN) are most valuable for the fainter and more extended 



TABLE 1 

Angular resolution for different radio observations. 



Instrument 


Wavelength 


Resolution 


VLBA 


43 GHz 


(0.21 X 0.43) mas 


VLBA (MOJAVE) 


15 GHz 


(0.6 x 1.3) mas 


VLBA 


2.3 GHz 


(7.5 X 3.9) mas 


VLBA 


1.7 GHz 


(8.0x3.4) mas 


EVN 


5 GHz 


(1.0x2.0) mas 


VLA 


15 GHz 


~ 0.13" x 0.12" 


VLA 


22 GHz 


~ 0.10" x 0.09" 


Note. — The beam 


sizes are given as full width at 



half maximum (FWHM). In case different resolutions 
have been used in the analysis the highest one is given 
in the table. 



HST-1 feature. 



VLBA 43 GHz — The 43 GHz VLBA (Napie r etal.|1994) data 
from 2006 through 2008 were collected as part of an effort 
to study the dynami cs of the inner jet near the launch region 
( [Walker et aLl[2008) l. The 2009 and 2010 observations were 
part of a project to fi nd and study anoth er VHE/radio flare 
like that seen in 2008 ( Acciar i et al.|2009) . The observations 
were made on the 10-antenna VLBA using a total bandwidth 
of 64 MHz. The data were reduced in AIPS following the 
usual procedures for VLBI data reduction including correc- 
tion for instrumental offsets using the autocorrelations, band- 
pass calibration based on strong calibrator observations, and 
correction for atmospheric opacity based on the system tem- 
perature data. The a priori amplitude calibration depended 
on the gains provided by VLBA operations, which are based 
on results from regular single-dish pointing observations of 
Jupiter, Saturn, and Venus averaged over many months. The 
images are based on data that are both amplitude and phase 
self-calibrated. The flux scale for each epoch was set by nor- 
malizing the self-calibration gain adjustments to the a priori 
gains for observations above 30° elevation on those antennas 
with good weather and instrumental conditions for that epoch. 
The flux densities typically accurate to within about 5%. 

Three VLBA flux densities are provided. The first is 
the peak brightness on the core in an image made with a 
0.21 x 0.43 mas beam. The second is the integrated flux den- 
sity within 1 .2 mas of the core and represents the total emis- 
sion from a region within a projected distance of 0. 1 pc from 
the presumed position of the black hole at the radio core 
(0.1 pc = 340 R s for a 3.0 x 10 9 M G black hole). The third 
is the integrated flux density from the jet over the region 1.2 
to 5.3 mas from the core. 

VLBA 15 GHz (MOJAVE)— VLBA 15 GHz observations from 



the MOJAVE program (Lister et al. 2009 1 have been analyzed 
to obtain core fluxes over 16 epochs from 2001.0 - 2011.0. 
The calibrated (m, v) data were retrievec 114 and re-imaged uni- 
formly with the final maps restored using a 0.6 ma s x 1.3 ma s 
beam (position angle = -11°) following Abdo et al. (20091. 
The peak core fluxes measured from the resultant maps span 
typical values of ~1.0— 1.2 Jy/beam (Fig. fTh with two notable 
peaks of ~1.5 Jy/beam recorded in early-2008 and mid-2009. 

VLBA 2.3 GHz — Observations at 2.3 GHz were made on 8 
and 18 April 2010 using 10 VLBA stations as a part of ex- 
periments BH163. Each session has a total on-source time of 
~15 minutes with the total bandwidth being 64 MHz. In order 

1 14 http://www.physics.purdue.edu/MOJAVE/ 
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to obtain better u, v coverage, the short scan blocks (~ 2 min- 
utes per block) were distributed uniformly at several hour an- 
gles. The initial data calibration was performed in AIPS based 
on the st anda rd VLBI data reduction procedures (the AIPS 
The amplitude calibration with the opacity cor- 



TABLE 2 

Parameters of the fit to the 2010 VHE flare 
data 



cookbook 115 



rections was applied using the measured system noise tem- 
peratures and the elevation-gain curves of each antenna.The 
data were fringe-fitted and then averaged at short intervals 
(every 5 seconds and 1 MHz in time and frequency domains) 
before the imaging process in order to avoid the smearing ef- 
fects due to the time and bandwidth averaging at the HST-1 
region. Th e images were made in DIFMAP software ( |Shep-| 
herd|1997 l with the iterative phase/amplitude self-calibration 



processes. The off-source rms noises in the resultant images 
are 0.4-0.5 mJy/beam. The peak flux densities are provided 
for the core region with the synthesized beam of 7.5 x 3.9 mas 
at -5°, and the integrated flux densities are provided for the 
HST-1 region. The errors in flux densities are assumed to be 
5% based on the typical VLBA calibration accuracy at this 
frequency. 

VLBA 1.7GHz— 19 epochs of VLBA 1.7 GHz observa- 
tions from -2005.0 - 2008.0 have been obtained in an ef- 
fort to monitor the evolution of the HST-1 knot following 
its bright ening at X-ray, optical, and radio frequencies ( Har- 
ris et al. 2009). The results of the first 9 of these o bserva- 
tions (progra ms BH126, BH135) were presented in [Cheung] 
et al. (2007), where apparent superluminal radio features in 
the HST-1 complex were discovered. For the additional 10 
observations (programs BC167, BH151), the identical cali- 
bration and imaging procedure was followed for this analysis. 
From maps restored with the same uniformly weighted beam 
(8.0 mas x 3.4 mas elongated north-south), we measured the 
peak core brightness (5% errors are assumed). Integrated flux 
densities (10% errors are assumed) for the HST-1 were mea- 
sured from naturally weighted images (11.5 mas x 5.5 mas) 
using a 80 mas x 50 mas box that covers the entire radio com - 
plex resolved in the VLBA images (cf., Cheung et al.|2007| l. 



EVN— M87 was observed with the European VLBI Net- 
work (EVN) at 5 GHz seven times between November 2009 
and June 2010 as part of a project aimed at correlating the 
radio and high energy behavior of the source. Given the inter- 
esting episodes of activity of the source during the campaign, 
a few observations were scheduled as ToO and do not have 
the same array configuration as the other ones. In general, 
6 to 11 telescopes participated in the observations, with base- 
lines ranging between less than 100 km (as provided by MER- 
LIN stations) and 7 000-9 000 km (as provided by the Arecibo 
and Shanghai stations, respectively). The observations were 
carried out at 5.013 GHz, divided in 8 sub-bands separated 
by 16 MHz each, for an aggregate bit rate of 1 Gbps. The 
data were correlated in real time at JIVE using the so-called 
e-VLBI technique, which provides a fast turn-around of the 
results (as w as the case near the Feb. 2010 event; see |Giro-1 
letti et al.|2010 1. Automated data flagging and initial ampli 



tude and phase calibration were also carried out at JIVE us- 
ing dedicated pipeline scripts. Data were finally averaged in 
frequency within each IF, but individual IFs were kept sep- 
arate to avoid bandwidth smearing. Similarly, the data were 
time-averaged only to 8 s, in order to avoid time smearing. 
Final images were produced in DIFMAP after several cycles 

115 http://www.aips.nrao.edu/cook.html 



Parameter 


Value 


Unit 


Fit range 


55290 - 55299 


MJD 




1.69 ± 0.30 


d 


tfecay 

r d 


0.611 ±0.080 


d 


to 


55295.954 ± 0.094 


MJD 


<D 


(2.01 ±0.15) x 10" 11 


-2 -1 
cm s 



Note. — The fit results in a ^ 2 /d.o.f. = 
10.02/11 with chance probability P = 0.53. The 
parameters are defined in the text. 

of phase and amplitude self-calibration. Various weighting 
schemes were applied to the data to improve resolution in the 
core region and enhance the fainter emission in the HST-1 
region. For the core, uniform weights were used obtaining a 
typical HPBW of 1.0x2.0 mas (in PA -20°) and providing the 
peak brightness. For HST-1, natural weights were used result- 
ing in a 7.5 x 8.5 mas HPBW; a resolved structure is clearly 
detected and integrated flux densities f or the full region ar e 
given. For additional details, please see Giroletti et al. (201 1 1. 



VLA — Data from the VLA archive for 10 epochs have 
been used, selecting observations performed in A-array at 15 
and/or 22 GHz. The typical angular resolution of the VLA 
is - 0.13" x 0.12" and ~ 0.10" x 0.09" at the two frequen- 
cies, respectively, which permits one to resolve the core and 
the HST-1 feature. Data reduction was carried out in AIPS 
in the standard manner: the flux density scale was tied to the 
main gain calibrator 3C 286 (which was used in all the ob- 
servations) using SETJY with the available model. Phases 
were calibrated using the compact source 1224+035, obtain- 
ing good solutions. A few final iterations of the phase and 
amplitude self-calibration have been carried out to improve 
the image quality. Integrated flux densities for the core and 
HST-1 have been derived using JMFIT. 

Add itional 15 GHz VLA data for HST-1 from Harris etaLl 
(2009) are also shown in the light curve. 



3. THE 2010 VHE CAMPAIGN 

The 2010 flare — During the 2010 VHE monitoring campaign 
two episodes of increased VHE activity have been reported 



(Mariotti & MAGIC Collaboratioh1|2010b| |Ong & Mariotti 
|2010| l: The first episode took place in Feb. 2010 where a sin- 
gle night of increased activity was detected by MAGIC (Fig. [2] 
bottom panel around -35 d; detection significance > 5 s.d.^ 
Follow-up observations by H.E.S.S. and VERITAS did not 
reveal further activity at VHE. The second episode took place 
in Apr. 2010 and showed a pronounced VHE flare detected by 
several instruments triggering further MWL observations. In 
the following, the discussion will concentrate on this second 
flaring episode. 

The VHE activity of this second flaring episode is concen- 
trated in a single observation period between MJD 55290 and 
MJD 55305 (-15 days; see Fig. [2] bottom panel & Fig.j3j. 
This time-period is exceptionally well covered with 21 point- 
ings by different VHE instruments, resulting in an observation 
almost every night. It should be noted that during nights with 
(quasi) simultaneous observations by different instruments, 
the measured fluxes are found to be in very good agreement. 

The detected flare displays a smooth rise and decay in flux 
with a peak around MJD 55296 (April 9-10, 2010). A peak 
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Fig. 2. — VHE light curve of M 87 of the flaring episodes in 2005 (top), 2008 (middle), and 2010 (bottom). Integral fluxes are given above an energy of 350 GeV. 
The lengths of the gray bars correspond to the length of the gray shaded areas in Fig.[T] A time of days corresponds to MJD 53460, MJD 54500, and MJD 55270 
for 2005, 2008, and 2010, respectively. Flux error bars denote the 1 s.d. statistical error. Horizontal error bars denote the time span the flux has been averaged 
over. Note that in the case of time spans longer than one night the coverage is not continuous. 



flux of ~ 2.5 x 10 _u cm _2 s _1 is reached, which is about a 
factor 10 above the quiescent flux level of the source. The 
data on the rising part of the flare indicates a steady rise. On 
the decaying side the situation is more complex: two nights 
after the detection of the maximum flux (MJD 55298) all three 
instruments measured a low flux compatible with zero, while 
in the following 3 nights (MJD 55299/55300/55301) a higher 
flux is detected by VERITAS. 

To derive the timescales of the flare the VHE light curve 
from MJD 55290 to MJD 55299 is fitted with a two-sided 
exponential function 



O = O n X e - ] '-' MAT with 



At = Ar nse for t < t 
At = Ar icc ^ for t > t 



(1) 



resulting in an excellent description of the data with a 
^- 2 /d.o.f. = 10.02/11 and a corresponding chance probabil- 
ity of P — 0.53. The resulting fit parameters are summarized 
in Tab. [2] Very similar results are obtained when (i) fitting 
a longer time span (e.g. MJD 55275 - 55305) or when (ii) 
adding a constant offset as a free parameter to account for 



a possible baseline flux, though the chance probabilities of 
these fits are slightly worse then for the original fit described 
above. 

Flux doubling times of T d = ln{2) x At of t™ 6 = (1.69 + 

0.30) days for the rising and T* cay = (0.61 1 ± 0.080) days for 
the decaying part of the flare are derived, which signifcantly 
differ by a factor 2.77 + 0.62. The variability timescales de- 
rived for the 2010 VHE flare are the most precisely measured 
VHE variability timescales determined for M 87 to date. Pre- 
viously detected VHE flares only allowed for rough estimates 
due to the variability pattern, sampling, and statistics. 

Comparison with the 2005 and 2008 flares — The VHE light 
curve from 2005, 2008, and 2010 around the flaring episodes 
is displayed in Fig. [2] In all three flares, similar flux levels are 
reached. The apparent timescales of the order of days are also 
very similar. The time period over which activity is detected, 
is comparable in 2008 and 2010 but is slightly longer in 2005, 
though the gaps in the sampling make any definitive statement 
difficult. During the 2010 activity period a pronounced flare is 
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Fig. 3. — VHE light curve of M 87 zoomed on the 2010 flare. Also shown are 
the results of the fit of an exponential function to the data. Error bars denote 
the 1 s.d. statistical error. The results of the fit to the data are summarized in 
Tab. [2] Further details can be found in the text. 

detected which is well described by an exponential behavior 
(see previous paragraph). While the 2005 flare is compatible 
with such behavior (the sampling is considerably worse than 
in 2008 and 2010), the 2008 activity state seems more erratic 
with several maxima over a similar time period as in 2010. 
Taking the best fit function to the 2010 flare as a template, 
large parts of the 2008 flaring activity are not compatible with 
the 2010 behavior (see the Appendix for details). 

VHE flare duty cycle — Given the typical length of a VHE 
flare of order a day, the nightly VHE y-ray flux measurements 
presented in this paper can be used to estimate the duty cycle 
of M 87 for VHE flares. For example , following the approach 
presented in |Jorstad et al. (2001 116 one derives a duty cycle 
of ~28%. However, due to the observing strategy followed for 
M 87 (observations have been intensified after the detection 
of a high state), the data-set is biased. The derived duty cycle 
is, therefore, likely overestimated and should be considered 
an upper limit. To enable a more unbiased view and to give 
an estimate of the possible uncertainties the duty cycle for 
a range of threshold fluxes is presented, calculated from the 
number of data points above the threshold flux divided by the 
total number. For a threshold flux of ^thresh = (0.5/0.8/1) x 
10~" cm _2 s _1 the respective duty cycles are ~14% / 7 % / 
4%. 

4. OPTICAL POLARIMETRY 

Radio and optical emission from regions with in the M 87 
jet show a high degree of polarization (Perr man et al.| 1999] >, 
which is indicative of non-thermal emission processes in 
highly ordered magnetic fields (synchrotron radiation). In 
general, polarization and especially changes in the polariza- 
tion can be utilized to investigate the magnetic field struc- 
ture of the emission sites and put additional constraints on 
the emission models. Recently, changes in the polarization 
(amplitude and angle) have been used to investigate the origin 
of the y-ray emission in several blazars ( Marscher et al.|2008| 
[20T0l[Abdo"et al.|2010b||Barres de Almeida et al.|201o| . The 

16 Calculating the number of observations with a flux a factor 1.5 above 
the error weighted mean (for this data set ~ 1.8 X 10 -12 cm" 2 s"'). 



observed pattern - a change in the HE/VHE flux coinciding 
with a rotation in the electric vector position angle (EVPA) 
and/or a change in the relative polarization - could result 
from: a non-axisymmetric magnetic field through which the 
emission region propagates, a swing of the jet across the line 
of sight, or a curved trajectory of the emission reg ion e.g. fol- 
lowing helical magnetic field lines or a bent jet ( Abdo et al. 
|2010bp . 



HST 2001-2008 — Regular optical polarimetry observations 
of the M 87 j et ha ve been performed with the HST from 2004 
to 2006 (Sec. 2.4 Perlma rTet al.|201 1) . The measured degree 
of polarization and the polarization angle for the core region 
and HST-1 are shown in Fig.|4j left panel. While both compo- 
nents show variable polarization their overall behavior is very 
different: during the flaring period (2004-2007) HST-1 shows 
a clear correlation of the fractional polarization with the total 
flux, with the fractional polarization ranging from 20 to 40 %. 
The EVPA remains almost constant ~ 7 - 8° off the direction 
of the jet (-69.5°; gray line in Fig. |4] lower panel). The nu- 
cleus, on the other hand, displays a highly variable EVPA with 
a lower overall polarization of 2 to 10 %. The EVPA changes 
are not correlated with the total flux or the fractional polariza- 
tion and appear, within the given sampling, erratic. The HST 
polarimetry data covers the 2005 VHE flare but, apart from 
the previously mentioned correlation of the HST-1 photomet- 
ric flux and fractional polarization, no correlations connected 
to the 2005 VHE flaring episode are apparent from the data. 

LT2010— Triggered by the detection of the 2010 VHE flare 
the LT started to take regular optical polarimetry observations 
of M87 in April 2010. The resolution of the instrument is 
not sufficient to separate the core and HST-1 and, therefore, 
one value for the combined core + HST-1 region is given. 
The LT detects clear variability of the degree of polarization 
of the core + HST-1 complex and marginal evidence (~ 2cr) 
for variability of the EVPA (Fig. |4j right panel). The polar- 
ization degree changes between ~ 2 and ~ 12% but without 
any apparent correlation with the VHE y-ray flux. The mea- 
sured EVPA indicates changes at the time of the VHE flaring 
episode from ~ -130° to ~ -60° (~ 70° in 3 days) and later 

settles to a stable EVPA of 20° about 50 days after the 

flare, though the only marginal indication for EVPA variabil- 
ity and the lack of continuous monitoring before and directly 
after the VHE flare makes it difficult to interpret these changes 
in terms of a continuous rotation of the EVPA. The observed 
polarization variability pattern is compatible with the pattern 
previously observed by HST for the core region and different 
to what has been observed from HST-1. This indicates that 
the 2010 LT data is possibly dominated by emission from the 
core. 

5. VHE FLARES & MWL CORRELATIONS 

Over the last 10 years M 87 has been extensively monitored 
all across the electro-magnetic spectrum from radio to VHE 
(Fig. [TV This large data-set can be used to investigate MWL 
correlations and thereby probe the origin of the VHE y-ray 
emission. In principle, many different physics processes could 
contribute to the production of VHE y-ray emissi on in M 87, 
e.g. annihilation of massive dark matter particles (Baltz et al. 
2000 ), cosmic ray interac tions in the extended radio lobes and 
the surrounding cluster (Pfrommer & EnBlin 2003), particle 
acceleration in the relativistic jets, etc. The detected short- 
term variability with timescales of the order of days and the 
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Fig. 4. — Optical polarization for different components of the M 87 jet as measured by HST and LT. The LT p olari zation measurement of the core + HST- 1 
complex is scaled up by a factor of 10 to compensate for the host galaxy emission in the measurements (see Sect. |2.4| for details). The horizontal gray solid line 
in the lower panel marks the jet direction (-69.5°) while the horizontal gray dashed lines are located at -69.5° ± 90" (180° distance). 



limits on the location place stro ng constraints on possible sce- 
narios (Aharoni an et al.|2006a| l, leaving the inner jet and the 
close vicinity of the SMBH as the most probable emission 
sites. In the following, the MWL behavior of the two most 
prominent features in the innermost structure of M 87, namely 
the HST-1 knot and the core, are discussed in the light of the 
VHE flaring activity. 

HST-1— Between 2001 and 2008 the first bright feature in 
the jet resolved by the HST in the optical, HST-1, unde rwent 
a spectacular flare detected in radio , optical, and X-rays ( |Perl- 
man et al.|2003| [Harris et al.|20 03). The flare displayed a rel- 
atively smooth rise over several years with a flux increase by 
more than a factor 50 in X -rays and optical. The flux peaked 
in the beginning of 2005 ( [Harris et aLl[2006l [2009] l around 
the same time when the enhanced activity level a nd the first 
short term v ariability had been detected at VHE (Aharonian 
et al. 2006a I. HST-1 has been discuss ed as a p ossible site 
for the VHE y-ray emission (e.g. |Stawarz et aL"1 |2006| |Oie 



ung et al |[2007l |Harris et al.||2009| l. While the size of H5T- 
1 as a whole is too large to account for the short-term vari- 
ability detected at VHE (following causality arguments), high 
resolution VLBA radio observations resolve HS T-1 into sev - 
eral, partially unresolved sub-structures ( |Cheung et al.|2007] l. 
These sub-structures also display apparent superluminal mo- 
tion up to 4.3 c + 0.7 c (see also Giovannini et al. 201TJ. In 
combination with the detected synchrotron X-ray emission 
and strong polarization of the radio-to-optical continuum, this 
indicates that efficient in-situ acceleration of the radiating par- 
ticles is taking place in compact sub-volumes of the HST-1 re- 
gion, characterized by well-organized magnetic field and rel- 
ativistic bulk velocities. On the other hand, during the 2008 
and 2010 VHE flares HST-1 was in a low flux state without 
pronounced activity at radio or X-ray wavelengths, thus disfa- 



voring it as the origin of the VHE y-ray emission during these 
episodes. 

Core — The direct vicinity of the SMBH and the jet base 
have been proposed as possible production sites of the VHE 
y-ray emission (e.g.|Reimer et al.|2004|[ Ghisellini et al.|2005 
Tavecchio & Ghisellini 2008, Neronov & Aharonian||2007 



Rieger & Aharonian 2008 ; Lenai n et al.||20 08 ; Bark ov et 
20101 IGiannios et al'|2010l|Levinson & Rieger 2011*03 



al. 
M~8T 

10 41 erg/s; Perlman et al. 



is only a weak IR source (vL v 
2001b) and, therefore, VHE y-rays are most likely able to es 
cape even from the close vicinity of the SMBH without suf- 
fering strong absorption due to yy-int eractions (Neronov & 
Aharonian 2007; Brodatz ki et al.|2011|) , although the debated 



origin of the observed IR photons and the poorly known struc 
ture of the accretion disk in the M 87 core, which both play a 
crucial rol e in this context, shou ld be kept in mind (see the dis- 
cussion in [^eung]eFin20(^ and jLrS~aLj2009). TheM87 
jet base has been imaged with VLBI with sub-mas resolution 
(see e.g. Ly et al.||2007 and references therein). It shows a 
resolved, edge brightened structure to within 0.5 mas of the 
core (0.04 pc) and indications for a weak counter-jet feature, 
suggesting that the SMBH lies within the central beam of the 
VLBI observations. At even shorter distances to the core the 
jet has a wider open ing angle, whic h is interpreted as the jet 
collimation zone ( |Junor et al.| 1999| >. 

In 2008, densely sampled 43 GHz radio observation of the 
innermost jet region s revealed a flare o f the radio core (flux 
increase of ~ 30 %; Acc iari et al.|20 09). At the same time, a 
flare at VHE and a subsequent enhanced X-ray flux from the 
core region followed by a sharp decrease were detected. The 
observed MWL variability pattern supported the interpreta- 
tion that the VHE y-ray emission likely origina tes from the 
close vicinity of the SMBH near the jet base (Acciari et al. 
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2009 ). The observed MWL behavior (VHE and radio flux, 
radio map) is well described by a simple, phenomenological 
model where the VHE fla res are associated w ith the injection 
of plasma at the jet base (Acciari et al. 2009 supporting on- 
line material). As the injected plasma blobs travel down the 
jet, they expand and become transparent for radio emission 
leading to the observed radio feature. 

In contrast, radio observations taken in 2010 contempora- 
neous with the VHE flare show no enhanced radio flux from 
the core region (Fig. [TJ: VLBA 43 GHz ToO observations 
triggered by the detection of the VHE flare indicate a sta- 
ble flux state of the core (1.2mas) and the inner jet (1.2 to 
5.3 mas) at the previously detected flux levels. In addition, 
EVN 5 GHz, VLBA 2.3 GHz, and MOJAVE 15 GHz measure- 
ments also show no indication of an enhanced radio flux state 
from the core in 2010. A direct comparison of these data with 
the 43 GHz data is difficult given the differenc e in resolution 

In the op- 
and HST-1 



and the missing overlap during the 2008 flare 
tical, only a combined measurement for the core 
with limited sensitivity by the LT is available, which does not 
indicate any strong activity in the two components in 2010. 

On the other hand, Chandra X-ray observations of the core 
show an enhanced flux ~ 3 days after the peak of the VHE 
y-ray emission in 2010 (see Fig. [TJ. The flux is enhanced by a 
factor ~ 2 for a single measurement and than drops back to a 
lower state less than two days later. The observed variability 
timescale is significantly shorter (by a factor ~ 10) then the 
shortest X-ray variability measured previously from the M 87 
core (< 20 days; Harris et al. 2009). Further details on the 
Chandra X-ray data from 2010 will be reported in a compan- 
ion paper (Har ris et al.|201 1| . 

It should be noted that the X-ray fluxes measured from the 
core during the time of the VHE flares in 2008 and 2010 are 
the two highest measurements since the start of the Chandra 
observations in 2002. This coincidence can be interpreted 
as indication that the VHE y-ray emission in 2008 and 2010 
originates from the X-ray core region. During the 2005 VHE 
flaring episode no enhanced X-ray emission from the core 
was detected. At that time, HST-1 was more than a factor 30 
brighter than the core region in X-rays leading to uncertain- 
ties in the flux estimation of the core (e.g. ' eat-thy-neighbor' 
effect; see Appendix A of Harris et al. 2009 [ ) which could sup- 
press the detection of a core flare. 

6. DISCUSSION 

General considerations — The VHE and MWL data in this 
paper can be interpreted in two fundamentally different ways: 
(1) each of the VHE flares detected originates from a different 
emission region and/or process, or (2) they have a common 
origin. Support for the former interpretation comes from the 
difference in the VHE light curves (variability pattern and the 
overall duration of each episode), as well as from the differ- 
ence in the apparent MWL correlations in the different years. 
In this interpretation, all previously derived models remain 
possible and an additional explanation for the 2010 flare has 
to be found. Observational support for the latter interpreta- 
tion, involving a common origin of the VHE flares, follows 
from the similarities in the detected peak fluxes, flux doubling 
timescales, and also in the spectral shapes of t he different 
VHE flares (for details on the VHE spectrum see Aharonian 
|et al.|2006aHAlbert et al.|2008aHAliu et al.|2011) . In this in- 

117 Noteworthy, one of the highest flux states in the MOJAVE 15 GHz data 
is measured in 2008 shortly after the VLBA 43 GHz observations ended. 



terpretation the only remaining MWL signature for the VHE 
flares is an X-ray flare of the core region (see previous section) 
locating the VHE y-ray emission site in the central r esol ution 
element of the Chandra observations (0.6 " ~ 50 pc) 



us 



In addition, it is not known whether the quiescent and flar- 
ing VHE y-ray emission are produced in the same region and 
by the same process. While, in principle, they can be pro- 
duced in two distinct locations by two different processes, the 
similarities of the VHE y-ray spectrum between the two states 
might suggest a common origin, or at least a very similar 
physical process involved. A common origin of the quiescent 
and flaring VHE fluxes is therefore anticipated below, even 
though different emission regions dominating the two states 
cannot be excluded. Interestingly, the HE spectrum measured 
by Ferm/-LAT joins smo othly with the V HE spectra derived 
for the quiescence state ( Abdo et al. 2009 ) possibly indicating 
a common origin and, therefore, HE- VHE flux correlations 
could be expected. Unfortunately, due to the low flux of M 87 
in the HE band and the limited sensitivity of ferm/-LAT, no 
conclusive statement on such correlation can be derived from 
the data presented in this paper. 

In the following, the discussion will mainly focus on the 
case where a common origin of all VHE flares in the X-ray 
core (0.6 " ~ 50 pc) associated with an X-ray flare is assumed. 

Characteristics of the VHE flares — The VHE flares are char- 
acterized by a variability timescale f var ^ 1 day, and a broad 
band power law spectrum extending from ~0. 1 TeV up to 
at least 10 TeV with photon index T ^ 2.3, weakly vari- 



2008a Aliu et 



able during flares (jAharonian et al. 2006a Albert et al. 

aLl|201 1|). The observed VHE flux during 



the flaring episodes reached similar maximum flux levels of 
®>o.35TeV - (1 - 3) x 10~' 1 ph cirT 2 s~ 1 . With the adopted 
distance of d — 16.7 Mpc and using the average photon in- 
dex T = 2.3, the corresponding isotropic VHE luminosity is 
Lvhe = 4nd 2 <D >£o E (F-l)/(r-2) - (0.8-2.4)x 10 42 erg s" 1 . 
If extrapolated down to the HE range (0. 1 Ge V) with the 
same photon index, the total y-ray luminosity of the flaring 
events would be even higher, namely Lhe-vhe ~ 10 43 erg s . 
This is a non-negligible amount given that the bolometric 
accretion luminosity of the M 87 nucleus is relatively low. 



10 42 ergs' 1 (jReynolds et al.||l996] |Di Matteo et al. 



|2003| |Kharb & Shastri 2004), and the total kinetic lu minos 
ity of the jet is also quit e modest, L\ ~ 10 44 ergs~' ( |Bick-| 



nell & Begelmann 1996; Owen et al j[2000| ). The efficiency 
of the VHE y-ray production in M 8 7 is therefore an issue, 
even taking into account order-of-magnitude dimmer quies- 
cence fluxes. More specifically, the observed VHE luminos- 
ity during the flaring events (which is approximately equal to 
the average/quiescence bolom etric y-ray lum inosity in both 
HE and VHE bands; see |Abdo et al.| 2009 ) is of the order 
of the accretion luminosity in the M 87 system, constituting 
at the same time about 1% of the jet total kinetic luminosity, 

ivHE ~ iacc ~ 0.01 XLj. 

Interestingly, assuming the observed emission is moder- 
ately beamed with a Doppler factor 5 of the order of a few 
and that, when viewed at smaller inclinations, the beaming of 
the nuclear jet in the M 87 system would be the same as the 
one deduced for blazar sources, namely 6* - 10-30, the 
isotropic VHE flaring luminosity of M 87 observed at smaller 
viewing angles would read L* F ^ (6* /6) 4 Lvhe - (10 44 - 



118 It should also be noted that direct observational limits on the extent of 
the VHE emitting region constrain its p rojected size to be < 14kpc centered 
on the radio core (Aharonian et al. 2006al. 
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10 46 )ergs _1 . This is in the range of VHE flaring luminosi- 
ties of blazars of the BL Lac type, for which low -power radio 
galaxi es like M 87 are believe d to constitute a parent popula- 
tion ( |Urry & Padovani| [T995 ). Moreover, in such a scenario 
the observed variability timescales are scaled down by the ra- 
tio of the Doppler factors At* = (6/6*) At ~ (0.5 - 0.2) At 
becoming of the order of a few hours, which is less than the 
timescale derived from the size of the Schwarzschild radii 
via causality arguments (~day). Variability timescales shorter 
then the causality timescales are also observed in flaring VHE 
blazars, where, in some cases, variability timescales down to 
a few minutes imply even higher Doppler factors of (9(100) 



(e.g. lAlbert et al.|2007| |Wagner|2008| [RTE S.S Collaboration 
|et al.|2010[ ). This agreement can be considered as support for 
a blazar-like origin of the y-ray emission in M 87. The VHE 
flux changes by a factor of 10 and a VHE flare duty cycle 
of (9(10%), as estimated in t his paper, w ould also be consis- 
tent with such a hypothesis ( Wagn erj2008| ). Note in this con- 
text that with relativistic beaming involved the power emitted 
during the VHE flares would be reduced with respect to the 
isotropic luminosity estimated above as L emi vhe - rj 2 L VHE , 
where Tj is the bulk Lorentz factor of the emitting region. 

Particularly interesting and constraining are the distinct 
timing properties of the one-day-long VHE flares of M 87, 
recognized clearly and quantified for the first time in this pa- 
per. These include repetitive outbursts (2005), erratic flux 
changes (2008), exponential flux increases/decays for well- 
defined isolated events, and significantly different rise and de- 
cay times (2010). If indeed the three flares share a common 
origin, every model dealing with the VHE y-ray emission of 
M87 has to be able to accommodate such a diversity in the 
observed variability patterns. 

Lastly, MWL correlations (or the apparent lack of such) 
have to be addressed as well, meaning in particular the emerg- 
ing connection between the VHE and X-ray bands with no ac- 
companying flux variations at radio or optical frequencies. It 
should be noted, however, that unlike the order-of-magnitude 
flux increases observed in the VHE regime, the observed vari- 
ability in the X-ray regime is of a much lower amplitude 
(flux changes by a factor of ~ 2 only), corresponding to a 
rather moderate X-ray core peak luminosity of the order of 
~ 10 41 erg s _1 . This value has to be taken with some caution 
though since, up to now, there are no truly simultaneous X-ray 
core observations during a VHE flare (i.e. during the night of 
the peak of the VHE y-ray emission). 

In the following paragraphs, existing models for VHE y- 
ray emission from the M87 core are briefly reviewed and 
discussed in the light of the new observational result pre- 
sented in this paper. The majority of models discussed have 
been published b efore the detection of M 87 by Fer/ra-LAT 
( |Abdo et ak1|2009| l. Models published before 2006 are also 
not constrained by th e VHE short-te rm variability discovered 
by H.E.S.S. ( jAharonian et al.|2006a) . 



Observations versus modeling — Let us start with the 'mis- 
aligned blazar'-type, or rather 'misaligned BL Lac'-type sce- 
narios (noting at the same time that the standard leptonic 'ho- 
mogeneous one-zone SSC approach has been often consid- 
ered to fail in explaining the observed VHE properties of the 
M 87 nucleus; but see also the discussion in Abdo et al. 2009 ). 
One of the fi rst models of this kind was discussed by |Reirrier] 
et al. (2004), and involved a mixture of hadronic and leptonic 



for the emission region size R ~ 10 16 cm, variability timescale 
f var ~ 1 day, and the expected moderate beaming, 6 ~a few. 
Designed to account for the non-simultaneous data known at 
that time, the HE/VHE continuum was dominated by the syn- 
chrotron radiation of protons and muons in strong magnetic 
fields. While proton synchrotron radiation, with a maximum 
intrinsic cut-off at ~ 0.3 TeV, appears unable to account for 
the measured VHE spectra detected up to ~ 10 TeV, syn- 
chrotron radiation from muons (produced in charged pion de- 
cays) and the pion cascade components could potentially ex- 
tend the VHE spectrum to higher energies. This would then 
imply pion-production losses to be at least of the same order 
of magnitude as proton synchrotron losses. The low energy, 
millimeter-to-UV part of the spectrum is ascribed to the syn- 
chrotron emission of primary electrons with the same injec- 
tion spectral index as the primary protons, and the high, non- 
simultaneously measured X-ray flux considered to be domi- 
nated by an additional component. A direct VHE/X-ray cor- 
relation with a high X-ray flux level is thus not expected in 

s uch a scenario. 

Georganopoulos et al. (2005) discussed a leptonic version 
or the blazar-type modeling of the M 87 nucleus, relaxing the 
assumption regarding th e homogeneity of th e emission re- 
gion. More precisely, Georganopo ulos et al.| considered the 
case of a relativistic jet decelerating substantially on sub-pc 
distances from the core, and showed that the velocity dif- 
ference between a faster and a slower portion of the outflow 
could lead to the enhanced inverse-Compton emission of the 
former one in the TeV range (due to a relativistic boost of the 
energy density of soft photons produced in one portion of the 
flow in the rest frame of the other), at a level allowing a fit 
to the VHE data. In the model, the core emission observed at 
longer wavelengths (including the GeV range) was predomi- 
nantly due to the slower and outer parts o f the decelerating jet. 
A related scenario was analyzed by Tavecchio & Ghisellini 
(2008), who assumed instead radial velocity stratification, in- 
volving a relativistic jet spine surrounded by a slower sheath. 



In the particular model fits presented by Tavecchio & Ghis 
|ellini| the observed radio-to-HE emission of the unresolved 
M 87 core was dominated by the radiative output of the jet 
spine, while the observed VHE y-ray emission was produced 
mainly in the jet boundary layer. Both models are charac- 
terized by compact sizes of the VHE y-ray emission region 
(R < 10 17 cm), moderate beaming, and no obvious — if any 
— correlation between the VHE and lower-frequency bands. 
Importantly, the VHE spectra calculated in the framework of 
both models were rather soft, due to the Klein-Nis hina and y 



ray opacity effects involved (see the discussion in |Tavecchio 



& Ghisellini|2008[). Multi-zone approaches, meaning a dou 



emission processes operating within its relativistic outflow at 
distances from the jet base small enough to assure R < c t vm -6 



bled number of the model free parameters, could possibly al- 
low for some modifications to the presented fits enabling one 
to accommodate some additional MWL constraints (for ex- 
ample the VHE/X-ray correlation), but the observed flat flar- 
ing spectra of M 87 in the VHE range will remain problematic 
in both cases. 

Blazar type models assuming a different structur e of the 
emission reg i on are discussed in Lenain et al. (2008]) and Gi- 
|annios et al. (2010). In those, the bulk of the observed emis- 
sion of the unresolved M 87 jet was proposed to be produced 
in extremely compact sub-volumes of the main outflow. In the 
scenario analyzed by |Lenain et al.| such compact sub-volumes 
were assumed to be multiple blobs ejected from the core with 
modest or highly relativistic bulk velocities into a cone with 
a large opening angle. In the framework of the model con- 
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sidered by Giannios et al. on the other hand, one is dealing 
with fast mini-jets moving in random directions within a rel- 
ativistic 'large-scale' outflow. Complex setups of the models 
resulted in the fact that different variability patterns and vari- 
ous variability timescales could, in principle, be expected, de- 
pending on a particular choice of the model free parameters 
(for which broad ranges of values could be considered). In the 
particular fits presented by Lenain et al. ( 2008| , the authors 
considered for example the radio-to-U V continuum of the 
M 87 core to be produced by the extended though still unre- 
solved portion of the jet, and argued that the synchrotron and 
inverse-Compton emission of tiny blobs moving within such 
a jet may account for the observed nuclear X-ray and y-ray 
fluxes, respectively, including not only the quiescent but also 
flaring states. The anticipated VHE/X-ray correlation with no 
accompanying flux variations at longer wavelengths is thus an 
interesting feature of the model. 

Yet the crucial assumption regarding linear sizes of the 
emitting blobs significantly smaller than the Schwarzschild 
radius of the M 87 black hole, R <& R§, might be regarded 
as questionable (see the discussion in Begelman et al. 2008 ). 
A possible physical justification for the formation of such 
ultra-compact sites of the enhanced energy dissipation close 
to SMBHs in AGN jet s, and M 87 in particular, was given by 
Giannios et al. ( 2010| >, who speculated that this is the mag- 
netic reconnection process which may trigger compact beams 
of plasma ('mini-jets'), which then propagate with relativis- 
tic bulk velocities within a strongly magnetized, extended 
and similarly relativisti c outflow. In their ap plication to the 
specific case of M 87, |Giannios et al. ( |2010| l demonstrated 
that the synchrotron and inverse-Compton emission of the 
reconnection-driven mini-jets could account for the observed 
X-ray and VHE n uclear fluxes, somewhat in analogy to the 
results obtained by |Lenain et al.| p008), again under particu- 
lar model assumptions regarding the model free parameters. 
In the end, not only the fast VHE variability and the flat VHE 
spectra, but also the VHE/X-ray correlation could be success- 
fully accommodated in the model. 

Possibly challenging for scenarios involving ultra-compact 
emission regions is also the required efficiency of the VHE y- 
ray production. As discussed above, the VHE flares observed 
in the M 87 system are characterized by an isotropic luminos- 
ity constituting at least 1% of the total kinetic power of the 
M 87 jet. Whether such a power can indeed be efficiently dis- 
sipated within tiny sub-volumes of the jet remains an open 
question, even taking into account the expected beaming cor- 
rections. 

Magnetospheric particle acceleration and emission models 
have also been invoke d to explain the observed VHE emis- 
sion in M87 (see e.g. Rieger |2011| for a review). Assumed 
to be operating close to the event horizon of the central su- 
permassive black hole, the anticipated variability timescale 
(light-crossing argument) could be as small t =s 2GMbh/c 3 - 
0.35 day, and thereby satisfy the observed variability con- 
straints. The observed radio- VHE correlation during the 2008 
VHE flare has been in terpreted to provid e additional support 
for such an approach ( Acciari et al.|2009) . Usually, efficient 
particl e acceleration i n these scenarios is related either to gap- 
type ( |Neronov~& Aharoni an|2007| |Levi nson & Riege rp011| l 
or centrifugal-type processes (Rieger & Aharonian 2008 ) oc- 
curring in the magnetosphere around a rotating black hole. 
While the former can be very efficient and lead to the onset of 
an electromagnetic pair cascade (triggered by the absorption 
of ambient photons up-scattered to high energies by electrons 



accelerating in the gap), the latter is less efficient such that the 
VHE spectrum is expected to be shaped by the ambient soft 
photon spectrum (see also Vincent & LeBohec 2010 1. Both 
approaches appear to be able to reproduce the observed hard 
VHE flaring characteristics, but may need an additional con- 
tribution to account for the HE continuum. 

Recently, a different type of modeling has been brought 
forwar d to explain the ob served VHE y-ray emission from 
M87: Barkov et al. ( 2010| l argued that interactions of a rel- 
ativistic and strongly magnetized outflow around the jet for- 
mation zone with a star partially tidally disrupted by the in- 
teraction with the SMBH can lead to a very efficient produc- 



tion of hadronic -originating VHE photons (see also Bednarek 



& Protheroe 1997 for a different version of the jet-star in- 
teraction model for blazar-type sources). The predicted ex- 
ponential character of the flux increase and decay during the 
resulting VHE event, a flat GeV to TeV y-ray spectrum from 
proton-proton interactions, as well as the accompanying X- 
ray flux enhancement due to the free-free emission of the 
shocked cloud of the stellar matter, are particularly interest- 
ing features of the model which should be kept in mind. 
Other scenarios could yet be considered in the same con- 
text as well, involving stochastic acceleration of high en- 
ergy particles within a turbulent accretion flow close to the 
ev ent horizon of a SMBH (in analogy to the model proposed 
by Liu et al. 2006 in the context of VHE y-ray emission 
of Sgr A*), or reconnection-driven impulsive acceleration of 
electron-positron pairs in a magnetized corona of the accre- 
tion disk, e. g., in analogy to the model proposed by Zdziarski 
|et al.| (2009) in the context of VHE observati ons of the G alac- 
tic X-ray binary system Cygnus X-l (see also de Gouveia Dal| 
|Pino et al.|20T0] >. 

HST-1 — The arguments against association of the detected 
VHE y-ray emission with the non-thermal activity of the HST- 
1 knot are twofold: First, rapid variability of the VHE con- 
tinuum involving day-long flares implies that the emission 
region size is significantly smaller than the inferred size of 
the knot. Second, apart from the 2005 event, no correla- 
tion between the VHE activity and the radio-to-X-ray syn- 
chrotron radiation of the HST-1 flaring region, variable on the 
timescales of weeks and months, has been established. The 
first argument has to be taken with caution, though. That is 
because the HST-1 flaring region is basically unresolved even 
for radio interferometers, as already noted before. More im- 
portantly, very recently it has been found that rapidly variable 
high-energy radiation can be generated far away from the cen- 
tral jet engine, within the outer parts of relativistic outflows. 
In particular, the analysis of the y-ray emission of a bright 
quasar PKS 1222+216 indicates that the day-long VHE flares 
in this object originate most likely in compact sub-volumes 
of a relativistic jet at parsec distances from the active cen- 
ter flAleksic et al.||2011||Tanaka et~aL|2011||Tavecchio et aT] 
201 l| l. Interestingly, analogous phenomena seem to occur on 



even larger scales as well. For example, Chandra monitoring 
of the Pictor A radio galaxy reveals that the synchrotron X-ray 
emission of the knots located at tens of kpc from the jet base 
is variable on the timescales of years. This is much shorter 
than the variability timescale expected following the causality 
arguments for a given jet radius at the position of the knot, 



which is of the order of thousands of years (Marshall et al. 



2010 1. Hence, it remains formally possible that also in the 



case of the HST-1 knot, located about lOOpc from the M87 
nucleus, rapid VHE flares are being generated on timescales 



16 



The H.E.S.S., MAGIC, & VERITAS Collaborations and the M 87 MWL Monitoring Team 



shorter than the ones expected for a homogeneous outflow. 

Yet, the indications found in this paper for correlation be- 
tween the VHE flares and an X-ray flux increases of the nu- 
cleus of M 87 provide observational support for the idea that 
the observed VHE y-ray emission, at least during the flaring 
episodes, is associated with the innermost parts of the jet or 
with the closest vicinities of the SMBH in the center of this 
radio galaxy. Again, in principal, it remains plausible that, 
while the unresolved core dominates the flaring states, the qui- 
escent VHE y-ray emission of M 87 is instead related to the 
HST-1 knot. However, above we have advocated for a com- 
mon origin of the quiescence and flaring VHE fluxes based 
on the spectral similarity of the two states. If the favored hy- 
pothesis is correct indeed, then interesting constraints on the 
physical parameters of HST-1 can be derived. That is because 
the inverse-Compton up-scattering of ambient photon fields 
to the VHE range by the high-energy electrons producing the 
observed radio-to-X-ray synchrotron radiation of the knot is 
at some level inevitable, with the efficiency depending most 
crucially on the unknown magnetic field intensity within the 
considered jet region. Since for the equipartition value of the 
jet magnetic field a relative ly intense VHE y-ray emission of 
HST- 1 should be expected (Stawar z et al.|2006 Cheung et al.' 



2007} , the upper limits for the radiative output of the knot in 



y-rays (following from the fact that the observed VHE flux is 
associated with the nucleus) implies a strong, possibly even 
dynamically relevant jet magnetic field at hundreds-of-pc dis- 
tances from the jet base. This should be then considered as an 
important support for the MHD models of AGN jets, like the 
one presented in the part icular cont ext of the M 87 core and 
HST-1 knot by (Nakamur a et al.[2010) , and discussed further 
from the observational perspective by |Perlman et al. (2011) 
and Chen et al. (201 1 ). Analogous constraints emerging from 
the VHE data ca n be investigated for the outer parts of the 
M87 jet as well (Stawarz et al. 2005; Hardcastle & Croston 
|20TT) H 1 " 

7. SUMMARY & CONCLUSIONS 

During a joint VHE monitoring campaign on the nearby 
radio galaxy M87 in 2010, a major flux outburst was de- 
tected, triggering further VHE and MWL observations. The 
coordinated observations led to the best-sampled VHE light 
curve during a flaring state from this source (21 observa- 
tions in 15 days), revealing a single, isolated outburst. The 
measured VHE light curve of the flare is well described by 
a two-sided exponential function with significantly different 
flux doubling times of r™ e = (1.69 + 0.30) days for the ris- 
ing and T * cay - (0.611 + 0.080) days for the decaying part, 
i.e. a difference of a factor 2.77 ± 0.62. This measurement 
provides the shortest and the most precisely determined VHE 
variability timescale of M 87 today, and the first detection of 
a significantly asymmetric VHE flare profile in the source. 

In comparison to previous VHE flares detected in 2005 and 
2008, the 2010 flare shows similar timescales and peak flux 
levels, but the overall variability pattern is somewhat different 
from the more extended periods of flaring activity with sev- 
eral flux maxima observed before, though the statistics and 
the sampling of the previous VHE flares limit a definitive con- 
clusion on that matter. From the VHE long-term light curve 
the duty cycle for VHE flares is estimated to be < 4 - 28%, 
depending on the assumed threshold flux defining a VHE high 
state. 

VLBA 43 GHz observations, triggered by the detection of 



the VHE flare, show no indications for an enhanced radio 
emission from the jet base in 2010. This is in contrast to ob- 
servations in 2008, where the detection of a radio outburst of 
the core contemporaneous with the VHE flare lead to the con- 
clusion that the VHE y-ray emissio n is likely prod uced in the 
direct vicinity of the SMBH (Acciari et al. 2009). Chandra 
X-ray observations, taken ~ 3 days after the peak of the VHE 
y-ray emission, show a high flux state of the core region in 
2010, supporting the interpretation that the VHE flare origi- 
nates from the innermost jet regions. 

The long term (2001-2010) light curve of M87, spanning 
from radio to VHE, is investigated for a common MWL signa- 
ture accompanying the three VHE flares. No unique signature 
is found. Observations of the jet component HST-1, which 
has also been proposed as a possible emission site of the VHE 
y-ray emission, show no enhanced activity in 2008 and 2010, 
disfavoring HST-1 as the origin of the VHE flares during those 
years. In 2008 and 2010 the VHE flares are accompanied by 
a high state of the X-ray core with a flux increase by a fac- 
tor ~ 2, while in 2005 the strong flux dominance (more than 
a factor 30) of the nearby X-ray feature HST-1 could have 
suppressed the detection of such an increase of the core emis- 
sion. Associating the VHE flares with the X-ray flares from 
the core places the emission site in the central resolution ele- 
ment of Chandra (0.6 " ~ 50 pc). 

Several models have been proposed to explain the observed 
VHE y-ray emission from M 87, most of which were based 
on the 'misaligned BL Lac' hypothesis. And, in fact, sev- 
eral observed properties of the source, including the broad- 
band y-ray spectrum, the multi-wavelength character of the 
VHE flares, and their overall energetics, could be consid- 
ered as support for the blazar-like models for the VHE y- 
ray emission. Yet, the particular one-zone or two-zone emis- 
sion models proposed involving only moderate beaming, both 
hadronic and leptonic, have difficulties in explaining VHE/X- 
ray correlated variability, for which indications are found in 
this paper, as well as the relatively flat VHE spectrum of 
M87 (photon index F ~ 2.3) extending up to ~ lOTeV. 
More complex blazar-type scenarios, involving ultra-compact 
emission regions and more substantial beaming, may, in- 
stead, face difficulties in accounting for the flares' energet- 
ics with the isotropic peak luminosity during the flares reach- 
ing Lvhe — (0.8 - 2.4) x 10 42 ergs~', which is of order of 
the accretion luminosity and about 1 % of the total jet power. 
Other non-blazar emission models discussed in the context 
of the VHE observations of M 87, which are quite success- 
ful in reproducing some of the established properties of the 
source, seem, on the other hand, challenged by the broad-band 
character of the y-ray continuum, as long as the GeV photon 
energy range constrained by the ,Ferm/-LAT observations is 
included. All in all, even though no emission model could 
be identified as the most plausible one, important limitations 
for most of them have been identified and novel observational 
constraints have been presented in this paper. Thus, the case 
of M 87 demonstrates the relevance of systematic, long-term 
and multi-wavelength monitoring of nearby radio galaxies in 
understanding the origin of the high-energy emission of radio- 
loud AGN. 
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APPENDIX 

QUANTIFYING THE DIFFERENCE BETWEEN THE 2008 AND THE 2010 VHE FLARES 

To investigate the compatibility of the 2008 and 2005 VHE data with the flare detected in 2010, the best fit exponential function 
to the 2010 data (Fig. [3] Tab.[2]i is taken as a template and compared with the 2008 and 2005 light curves. The position of the peak 
fo is varied (all other parameters fixed) and the resulting ,y 2 /d.o.f. and probabilities are calculated for the data within a time span 
of length of the original fit length: fg - 7 < t^ ata < to + 2. The results for the 2010, 2008, and 2005 data around the flaring episodes 
are shown in Fig. [5] with the gray band again denoting the observation period with the flaring episode in each year. The statistics 
in 2005 are low resulting in several good probabilities (high probabilities > 0.1) over the flaring episode. For most parts of the 
2008 flare, the probability for the 2010 fit describing the data is well below 1%. The only notable exception is the beginning 
of the episode around MJD 54496 where the first two data points are reasonably well described by the function resulting in a 
probability of 0(5%). Note that the peak flux from the fit (~ 2 x 10~" cm _2 s _1 ) is similar to the peak fluxes measured in 2008. 
Varying the fit parameters in their errors does not significantly change the results. For comparison, the same study is shown for 
the 2010 data where, at the best fit position from the fit, the probability is ~70%. 
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Fig. 5. — VHE light curve (top), reduced chi-square^ 2 /d.o.f. (degree of freedom) (middle), and probability (bottom) vs the position of the peak to of the best 
fit function to the 2010 VHE data for the 2005 (left), 2008 (middle), and 2010 VHE data (right). Values are calculated for data points inside a time-span of 
to - 7 < ((j ata < to + 2. The gray bands mark the same time period as in Fig.[T|and|2] denoting periods of increased VHE activity. Dashed vertical lines mark the 
positions of VHE flares. The horizontal line in the bottom panels marks a probability of 5%. 
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